UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTINE FA 1 - 122 



TITLE 

NUMBER OF WORDS 
TYPE 



TEMPORARY STORAGE 
DURATION 



ACCURACY 

PAEAMETERS 
LOCATION 



Second Order Linear Differential Equation with 

Two Point Boundary Conditions (DOI Only) 

101 

Closed auxiliary for code Al^ Code XAl, the floating 

point constant listing auxiliary, must also he used. 

The routine is entered as follows s 



p 


50 gF 
50 pF 


p + 1 


26 xF 



where g is thfe l«^cation of the auxiliary routine and 

X is the location of Code FA 1. 

See parameters , 

T = (n + l) (t + 400) + ^00 approximately where T is 

the total time in milliseconds, t is the time in 

milliseconds of the auxiliary routine and n is the 

number of sub^intermls into which the range of x 

Is diTided* 

9 decimal places or less, depending on the condition 

of the given equation* 

Preset parameters are S3> S4, S5, S6, ST« 



CONTEHT 
OOF 00 mF 

OOF OOjfe 
OOF 00 kF 



USE 

m is the location of the floating point 
accumulator . 

^^is the location of code Al 
Locations :k + *) (j =0, .•*, 2n +5) are 
used for temporary storage. The 
aolutlpns y. ( i = 0, ^.. , n) will 
appear in the n + 1 successiYe locations 
k:+ n + 5 to k + 2n + 3^ 



-2^ 



OOF 00 qF 



7 



ro 



OOF 00 rF 



EROGRAM PAmMSTEl 

8 OOF 00 aF 



L©catiGns q. + j (j = 0, 1,..*_^10) cgb-^ 
tajba the difference eqiiation coefficients 
a^> h^, e,> d.^ -a^e differential eq-uation 

coefficients, A. * B. , C, , D. * ( i = 0, 

i 1 1 1 

If *-»M , n), the increment A (x) (8s6), 
aM a temporary storage (9S6), n?he 
differential eqiiation coefficients. A., 
B^, 0^^ D., calculated by the atJXiliary 
rdutine^i must "be stored in locations ks6y 
5S6, 6s6, and 7S6 respectively » 
liOeations r+j (j = 0, 1^ *, *, 7) must 
contain the length of the range of x, L> 
and the 'boundary eqiaation coefficients 
E, F^ G, H, K, M e^cpressed as floating 
point numbers, at the start of code 
FA 1* The quantities Ly E, F^ G^ H^ K 
and M are to be stored in locations ST, 
ISJ^ 2S7, 3ST> ^S7, 3S7s and 6s7 re* 
spectiYely, 



DlSCSRIPnOl 



n is the number of sub*ihtervals into 
vhich the range of x is divided* 
Over an interval ©f length L at eq.ually spaced points, 
x^(i = 0, 1, ,•* yn)y this routine will give the 
solution of the differential equation. 



Ay" + By« -f Qr = D^ 
vith the given boundaiy conditions, 

E^y* + Fy = G]. 



x=x. 



and 



Hy^ + Ey =. M] 



x=x - n^+x- 

m « 



^^ 



■where A, B, C, D, E, F, G, H, K and M are funetions of x^ the ludependemt 
variable « 

M open subroutine^ written using InterpretiTe orders, 
Biust be provided starting at location g. It Is to evaluate the coefficients 
A, B, e and B for a given x^ x < x < n Ax + x , of the differential 
ecjuation. These calculated coefficients A, B, C and B must be stored 
successively starting at Location kUS* ^is auxiliary subroutine must be 
left mth the word 

8k IF 

83 (address of the 8th word of 0ode Eal) F* 
HO!EE: Ax is available for use in the auxilial^ subroutine, 

Ax is stored in 8s6» 

MEMOD: This routine first calculates the coefficients of the 

difference ecLuaticats, 

Vh ^ Va-l = \ 

from the coefficients of the differential equations. The following 
equations are used to calculate a.> b., c and d (O < 1 < n): 

a. = A. + B. Ax/2 
1 ^ 1 ^ 



2 

b. «= C. • Ax - 2 A. for i 5^ 0, 

1 1 1 '^ * 



n 



^ = A^ - B^ Ax/2 Ax = L/n 

2 



d. = D. Ax' 

1 1 



A- 



^0 = -^0^0 






a = 
n 



b = (KB - C^H ) A x72 + A H + A K A 



2 

X 



n n n n n n n n 



c = -A H 
n n m 



d = A M^ • A X + (B M - H D ) A x^/2 

At each step, these coefficients are stored in the four locations starting 
with location q. Then this routine evaliiates the following formulae; 



9 



1 11 *'i-l 



1 



d. - c. z: , /r. 

1 i 1-1 ' r 1 



nr; - a./p^ 

When i = 0, ..* ^n and 'V^^ = = T^ . p is calciilated first and stored 
in location k + 1. Then ll^ is calculated and stored in location k + n + 3 
Finally -K^ is calculated and stored over i^ in location k + 1. Diiring 
the next time thru the loop^^^ t and ^ are ealciilated and stored in 
locations k + 2^ k + n + 4 and k + 2 resprectively. This process is con- 
tinued until ^^, f and n^ have been calculated and stored. At this 
' ■" n n 

stage, the "^ are stored in locations k + 1 thru k + n + 1 and the Z", 

1 1^ 

are stored in locations k + n + 5 thru k + 2n + 3. When this has been 
done, the routine calculates the solution using the formula 



-5- 



■where y ^ =0. These solutions^ y. are stored over the '?:^ in locations 
k + n + 3 thru k + 2n + 5* 

UOTE: Incorporated ¥-ithin the routine is a test to determine 

if 0. =0. When B . =0^ l/4 x 10 is stored in place of this 

The method was adapted hy S. Gill from a process used 
"by L. Fox (L< Fox and H. H. Eobertson| Automatic Digital Computation 
Symposium, National Physical Lahoratory, England^ March 1955* Paper 19^ 
pp. 5;, 6). 



RTs 9/27/55 
DAT E Fehruiary 16, 1954 

PROGRAMMED BY E. Polivka 

APPIOVID BY^ 



■j;0mI^ 



RP: mge 
9/ST/55 



^ LOCATIOir 


OEDER 


BOns PAGE 1 FA^ 





15 29F 
42 97L 






Set lip link and address of 


1 


10 20F 
42 7L 






aiixiliary siibroutine 


2 


36 83L 
26 85L 






5 


41 S5 
50 3L 






k 


26 s4 








Bkss 








5 


8s as6 








6 


85 ST 

86 8s6 
Bs as6' 






-IJalcijIate A x 


7 


8k IF 
83 (g)F 


Ftqie 33 j 

34 
Br 1 




to a-uxillarj sti"br0iitiiie 


S 


te 4S6\ 


By 74, 78 




replaced hy 8k IF 


9 


\misij 

8S S6 
81 5S6 






83 (54 )L 
Store a^ 


XO 


8T 2ST 
8s 9S6 








11 


85 6s6 








12 


8T 1S7 
84 9S6 
87 8s6 








13 


87 8s6 






Caldiaate the ialtial difference 




87 m-f5+00 






14 


80 s6 
8s 9S6 






eqiaatioa coefficients > a^, b , c 
and d^ 


15 


85 4s6 
87 8s6 






u 


16 


87 2S7 










84 9S6 




_! 





LOCATION 


ORDER 


NOTES PAGE 2 


. 17 


8S 1S6 






store b 




8k OF 






18 


8s 2S6 






store c^ 





81 5S6 






19 


87 3S7 

8s 9S6 








20 


85 1S7 
87 7S6 








21 


Qk 9S6 
87 8s6 








.22 


87 8S6 
87 N+5+00 






- 


23 


8s 9S6 
85 hs6 








24 


87 8s6 
87 3S7 








25 


84 9S6 
8s 3S6 






store d 


26 


8 J jkL 








OK SB 




— 




27 


81 2S6 
07 S5 


From 64^53 






28 


84 1S6 
8j 69L 








29 


OS 1S5 






store ». 

Ca Icii inte ff. s t., and T/' * 


30 


81 2S6 

07 (k+H+2); 


F 






84 5S6 








31 


06 1S5 


/ 






32 


OS (k+n+3) 
85 S6 


? 




store ^T 




06 1S5 









FA 1 



LOCATION 


ORDER 


NOTES PAGE 3 


53 


OS 1S5 
03 7L 






storel/t 
1 


5^ 


8j 77L 








83 (7)L 


by 82 






35 


Vps S6/ 






replaced by 8K IF 
82 6kL 


36 


8i ks6 
87 4S7 








37 


8s 2S6 
81 6s6 








38 


87 i^-S7 
8s 9S6 








39 


85 5S6 
87 5S7 








i^o 


8k 9S6 
87 8s6 








in 


87 8s6 
87 H+5+00 








k2 


80 2S6 
8s 9S6 








^3 


85 hs6 
87 8s6 








kk 


87 5S7 
Qk 9S6 








h^ 


8S 1S6 
81 ii-S7 






Calculate to the difference 

equation coefficients a ^ b 

n n 


k6 


87 7S6 
8s 9S6 






c and d . 
n n 


kj 


85 5S6 
87 6S7 








kQ 


Qk 9S6 










87 8s6 









FA 1 



LOCATIOl 


QRBER 


HOTES PAGE h 


h9 


87 8s6 
87 S+5+00 






\. 


50 


8S 9S6 
85 i^s6 








51 


87 8s6 
87 6s7 








32 


&k 9S6 
8S 5S6 








53 


8k if 

83 27L 








^h 


85 5S6 
87 8s6 


from 8 






55 


87 IT+5^G0 
8S 9S6 








56 


m ks6 
8s s6 






store a. 


57 


81 9S6 








58 


8s 2S6 
85 6s6 






store 

calculate the difference 


59 


87 8s6 






-equation coefficieiits a. 5 h.y 


6g 


87 8S6 
8g ks6 
80 ij-s6 






c and d. i ^ O.n. 
i 1 ^ 


61 


8s 1S6 
85 7S6 






store "b. 
1 


62 


87 8s6 
87 8s6 




( 




63 


8s 3S6 






store d. 
1 


64 


8k IF 
83 27L 
IK S8 


From 35 


^S 




' 65 


11 (k+22i+3)F 








17 fi4-ii)F 1 







FA 1 



LOCATION 



OBDER 



lOIES PAGE 5 



FA 1 



66 

6? 
6B 

69 
70 

71 
72 

75 
jk 

75 
76 

77 
78 

79 
80 
81 



14 (k+2n+2)F 

15 (k+2H+2)F 
IL 1022F 

15 65L 
8 J 94l 
32 94l 
L3 S3 
52 7OL 
26 2%4 
19 IF 
40 S3 
L5 75L 
itO 1S3 
26 29S4 
00 F 
00 28F 
L5 76L 
40 8l 
26 29S4 
00 7L 

Sk rr 
83 54l 

^9 S3 
L5 93L 
% 8l 
LI 98L 
kO 98L 
56 29S4 
L5 lOOL 
kO 35L 
L5 93L 
42 34l 



from 28 



froffi 26 



From 34 



-Calculate the solutions, y 

i 



iJest Tor 



^' 



= 



Replace word 8 
iy 8k if 83 5^L 



S"wltclij 1st time go to part of 
code to a 
2nd time 5 



code to calculate a ^ "b > c « d . 

n n n^ n 



Go to part of code to calculate 
the J. 



LOCATION 



OEDER 



HOTES 



PAGE 6 FA 1 



82 

m 

85 
86 

8? 
88 

89 
90 
91 
92 

93 
94 

95 



9T 



26 2984 
00 F 
L5 5F 
1,4 8f 
42 65L 
F5 65L 
L4 TQL 
42 92L 
GO 20F 
k6 50L 
F5 65L 
f4 TOL 
h2 3IL 
L5 92L 
l4 8f 
42 66l 
00 20F 
46 66l 
F5 66l 
Op 20F 
46 65L 

41 (k+ii-j-2)F 
26 5L 
00 35L 

Id 98L 
4o 8l 
L5 99L 
4o 3% 
L5 T5L 

42 3^L 
52 971. 
22 ( )F 



From 2 



B^r 85 



By 



vaste 



Set addresses in routine dependent 



an tlie program parameters . 



Reset words 



8^ 3h, and 55 



link 



LOCATIOI 



ORDER 



NOTES 



PAGE 7 



FA 1 



98 

99 
100 



85 hs6 
87 1S7 
8k of 
8s b6 
8k IF 
82 6kL 



ETs 
SEPT* 28, 1955 
MGE 



